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The recently demonstrated control over light distribution through turbid media based on real¬ 
time three-dimensional optoacoustic feedback has offered promising prospects to interferometrically 
focus light within scattering objects. Nevertheless, the focusing capacity of the feedback-based 
approach is strongly conditioned by the number of effectively resolvable optical modes (speckles). 

In this letter, we experimentally tested the light intensity enhancement achieved with optoacoustic 
feedback measurements from different sizes of absorbing microparticles. The importance of the 
obtained results is discussed in the context of potential signal enhancement at deep locations within 
a scattering medium where the effective speckle sizes approach the minimum values dictated by 
optical diffraction. 


Light scattering in nanoscale heterogeneities is respon¬ 
sible for the optical opacity of materials, in particular 
biological tisues [H. Even though photons can gener¬ 
ally penetrate scattering tissues up to a depth of several 
centimeters Q], progressive changes in their propagation 
direction impede the ability to trace the exact trajectory 
of light and hence see through this type of samples. The 
light intensity distribution caused by photon interference 
is however not entirely random but determined by the 
type and location of optical scatterers, which results in a 
characteristic speckle pattern Q. By controlling the in¬ 
cident light wavefront, it becomes possible that photons 
positively interfere at specific locations, thus making fo¬ 
cusing through scattering objects theoretically feasible. 
Wavefront shaping techniques have recently emerged as 
a promising tool to image behind strongly scattering sam¬ 
ples, with a handful of potential applications foreseen [J- 
ISl]. Of particular importance is the possibility to image 
deep into biological tissues, where scattering limits the 
penetration of the modern optical microscopy approaches 
to a depth of a few hundred microns [1^. In the so- 
called near-infrared window, diffuse light can penetrate 
up to several centimeters into highly scattering tissues 0 , 
thus wavefront-shaping-based techniques can potentially 
enable high-resolution bio-optical imaging at centimeter- 
scale depths, ultimately limited by light absorption only. 

Wavefront shaping methods are based on the spatial 
modulation of the phase of the incoming light, which is 
commonly done by means of a spatial light modulator 
(SLM) device [ij. The optimal phase mask at the SLM 
can be determined from readings of the light intensity at 
the target point (s). The maximum signal enhancement 77 
at a given point (focus) depends on the spatial resolution 
of the feedback method via the following relation [llj | 


A^slm 

‘^^modes 


( 1 ) 


where A^modes is the number of optical modes generating 


the measured signal and Aslm is the number of degrees 
of freedom (pixels) of the SLM. The number of optical 
modes generally refers to the number of speckles enclosed 
in the area (or volume) that can be effectively resolved 
with the technique employed for measuring the light in¬ 
tensity. 


The feasibility to focus through scattering objects with 
wavefront shaping techniques was first demonstrated by 
using the two-dimensional image acquired with a CCD 
camera as a means to determine the light energy distribu¬ 
tion . Application of this approach to fluorescence mi¬ 
croscopy was further showcased Optical approaches 
are however not suitable for providing an efficient feed¬ 
back on the light distribution inside scattering objects as 
the spatial resolution quickly deteriorates in the diffuse 
propagation regime 0 - Thereby, alternative approaches 
based on a combination of optics and ultrasound have 
been recently suggested as a means to provide high reso¬ 
lution feedback from a scattering medium. For instance, 
the frequency (wavelength) of the light beam can be mod¬ 
ulated at a given location with focused ultrasound so that 
the phase of the modulated wavefront can be holograph¬ 
ically recorded after propagation through the object. In 
this case, by means of the phase conjugation, the light 
beam can be directed to a certain ultrasonically-tagged 
position 12 Alternatively, wavefront shaping based 
on the optoacoustic feedback may offer important ad¬ 
vantages to control the scattered light distribution. In 
this approach, the feedback signals can be in principle 
collected at a very high rate, ultimately only limited by 
the ultrasound propagation speed, and hence provide a 
very fast control mechanism over SLM mask optimiza¬ 
tion. Acquisition of signals with transducer arrays com¬ 
bined with graphics processing unit (GPU)-based recon¬ 
struction procedures further enable real time image ren¬ 
dering in two- and three dimensions 0 , which can even 
exceed the currently available refresh rates of standard 
SLMs 0 . Thereby, both single-element transducers and 
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transducer arrays have been suggested for optoacoustic 
wavefront shaping 17H21j. Moreover, the high optoa¬ 
coustic resolution enables a high light intensity enhance¬ 
ment according to Eq.[I] which can be further enhanced 
by conveniently selecting the detection frequency band 
or by using non linear effects [1^. Finally, light ab¬ 
sorption essentially occurs in any substance in nature, 
which makes optoacoustics arguably the most versatile 
feedback approach in terms if imaging contrast. 

Typically, the photothermal energy conversion lead¬ 
ing to optoacoustic signal generation takes place in light 
absorbing molecules. The size and average separation 
between the molecules is generally much smaller than 
the speckle dimensions. For example, organic molecules 
have typical lengths of a few nanometers, while the min¬ 
imal speckle size is determined by the optical diffraction 
limit, i.e. A/2, where A is the wavelength. For more 
than one absorber per speckle, A^modes in Fq. [prefers to 
the number of speckles in the region that can be resolved 
with the method employed to measure the light intensity. 
However, larger particles sparsely distributed across the 
region of interest can also be used to generate the op¬ 
toacoustic signals, in which case the number of speckles 
may be generally higher than the number of absorbers. 
In this case, A^modes in Fq.[T]is given by 


N 


^modes — A^sp.ii 


( 2 ) 


being N the number of particles enclosed in the resolved 
region and A^sp,i the number of optical speckles in each 
particle. The ratio between the sizes of the speckles and 
the particles becomes then a crucial factor in determining 
the potential enhancement achieved, which is also con¬ 
ditioned by the number of particles required to achieve 
detectable optoacoustic signal levels. Proper selection 
of the absorber characteristics is therefore of key impor¬ 
tance for achieving favorable wavefront shaping results 
with the optoacoustic feedback approach. 

A previously reported volumetric (three-dimensional) 
optoacoustic wavefront shaping platform [2lj was used 
herein to experimentally investigate the effects of the 
particle size on the light intensity enhancement. The 
lay-out of the set-up employed is depicted in Fig. 1. 
Basically, a phase-only SLM consisting of a liquid crys¬ 
tal on Silicon (LCoS) microdisplay with 8 fim pixel 
pitch (PLUTO-BB II, Holoeye Photonics AG) was used 
to control the light wavefront distribution after passing 
through a ground glass diffuser (Thorlabs DGI0-I20). 
The light enhancement was targeted at polyethylene mi¬ 
croparticles with approximate diameters 100 /rm (Co- 
spheric BKPMS 90-106) and 200 /rm (Cospheric BKPMS 
180-210) as well as Carbon microspheres with a diameter 
around 400 fim (SPI-Supplies) embedded into an agar 
phantom. The microspheres were excited with the out¬ 
put beam of a frequency-doubled Q-switch Nd:YAG laser 



FIG. 1. Lay-out of the experimental set-up. 



FIG. 2. (a-c) Photographs of the agar-embedded microparti¬ 
cles with diameters 100 /rm, 200 /rm and 400 /rm, respectively; 
(d-f) Optoacoustic images of the corresponding phantoms ob¬ 
tained with a constant phase value of the SLM pixels (initial 
iteration). Maximum intensity projections along the depth 
direction are shown; (g-i) Gorresponding optoacoustic images 
obtained after the phase of the SLM pixels was optimized to 
deliver maximum light fluence to the particles labelled ’A’, 
’D’ and ’G’ respectively. 


(Lab-190-30, Spectral Physics) operating at 15 pulses per 
second, which was collimated and horizontally polarized 
before being directed onto the SLM. The generated op¬ 
toacoustic signals were simultaneously collected with a 
spherical array of piezoelectric transducers 2^ , while the 
three-dimensional optoacoustic images representing the 
light absorption distribution in the phantom were recon¬ 
structed for each laser pulse using a graphics processing 
unit (GPU) [il. 
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The reconstructed image was used as a feedback mech¬ 
anism for a genetic algorithm that optimizes the phase 
mask at the pixels of the SLM based on the maximiza¬ 
tion of a defined cost function 21|. In the experiments, 
a mask of 20x20 phase values was optimized by grouping 
the SLM pixels accordingly. The speckle diameter was 
set to approximately 27 /rm by controlling the distance 
between the particles and the diffuser. The cost function 
was defined as the maximum image value in the volume of 
interest (VOI) enclosing a certain microsphere, such that 
the genetic algorithm converges to an SLM mask that 
focuses light at this particular absorber. The genetic al¬ 
gorithm performs the optimization as described in 0. 
Specifically, a population of 20 phase masks evolves with 
a refresh rate of 25% and a mutation rate that ranges 
between 10% and 1.25% of the pixels of the phase mask. 

Actual photographs of the phantoms are depicted in 
Figs. 2a-c for microspheres with diameters 100, 200 and 
400 ^m, respectively. The maximum intensity projec¬ 
tions (MIP) along the depth direction of the correspond¬ 
ing reconstructed optoacoustic images are subsequently 
shown in Figs. 2d-f. The MIPs in Figs. 2d-f are obtained 
by setting a constant phase value in the SLM, which rep¬ 
resents the initial iteration of the genetic algorithm. The 
corresponding MIPs of the images rendered after 2500 
iterations of the genetic algorithm are displayed in Figs. 
2g-i for cost functions defined to maximize the light in¬ 
tensity on particles labeled A, D and G, respectively. The 
evolution of the cost function (maximum signal intensity) 
as a function of the number of iterations in the genetic 
algorithm for all particles labeled in Fig. 2 is displayed in 
Figs. 3a-c. Fig. 3d shows the signal enhancement for all 
the imaged particles as a function of their diameter. The 
enhancement is estimated as the maximum value of the 
cost function for the given particle after normalizing it 
to the average value for the first 20 iterations, which cor¬ 
respond to the initial random population of the genetic 
algorithm. The small enhancement for the 400 /rm parti¬ 
cles can be attributed to signal fluctuations approaching 
the noise levels, so that it was not possible to focus light 
in those cases. As expected, the enhancement is maxi¬ 
mized for smaller particles containing a lower number of 
speckles. 

The results showcased indicate the importance of 
the absorber dimensions used to provide feedback in 
optoacoustics-based wavefront shaping, which ultimately 
determines the focusing (light enhancement) capacity 
of this approach. The maximum signal enhancement 
achieved is of particular relevance when eventually aim¬ 
ing at focusing light inside scattering samples, e.g. bio¬ 
logical tissues, where the speckles at depths beyond the 
transport mean free path are expected to have charac¬ 
teristic lengths approaching A/2. In order to generate 
efficient feedback from absorbers in the sub-micron size 
range, the optoacoustic imaging system would need to 
have an efficient detection bandwidth in the hundreds of 
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FIG. 3. Value of the cost function as a function of the number 
of iterations in the genetic algorithms for the particles with 
diameter 100 fj,m (a), 200 /xm (b) and 400 /rm (c). The light 
intensity enhancements achieved as a function of the particle 
diameter are displayed in (d). 


MHz range. However, the image artifacts and resolution 
degradation associated with acoustic heterogeneities and 
attenuation of the high frequency components in the im¬ 
aged medium are expected to impose hard limitations 
on the achievable spatial resolution when using the op¬ 
toacoustic imaging feedback. In biomedical optoacoustic 
imaging, the practically achievable resolution is g enerally 
in the range of 1/200 of the imaging depth 23 . Thus, 
the realistically achievable spatial resolution at a depth 
of 5mm is expected to be in the range of 25 /rm. As a 
result, approximately 250000 speckles will be accommo¬ 
dated in each resolution-limited voxel, for which no sig¬ 
nificant light intensity enhancement is feasible. Higher 
efficiency of the light enhancement can be potentially ac¬ 
complished when properly selecting the absorbing parti¬ 
cles and their concentration so that a relatively low num¬ 
ber of light speckles is enclosed within the resolved vol¬ 
ume while simultaneously attaining sufficient detection 
sensitivity. Yet, focusing light within biological tissues 
still remains a challenging topic due to additional fac¬ 
tors, such as the millisecond-level speckle decorrelation 
times. 


In conclusion, the size of light absorbing particles is 
designated to play an essential role in volumetric wave- 
front shaping using optoacoustic feedback. The experi¬ 
mental analysis performed in this work may help to pro¬ 
vide general guidelines for optimizing the particle size for 
achieving efficient light enhancement at locations exhibit¬ 
ing small speckle size, e.g. beyond the transport mean 
free path depth within turbid tissues. 
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